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Non-linear Gradient Algorithm for Parameter Estimation: Extended 

version 


Juan G. Rueda-Escobedo and Jaime A. Moreno 


Abstract —Gradient algorithms are classical in adap¬ 
tive control and parameter estimation. For instantaneous 
quadratic cost functions they lead to a linear time-varying 
dynamic system that converges exponentially under persis¬ 
tence of excitation conditions. In this paper we consider 
(instantaneous) non-quadratic cost functions, for which 
the gradient algorithm leads to non-linear (and non Lip- 
schitz) time-varying dynamics, which are homogeneous in 
the state. We show that under persistence of excitation 
conditions they also converge globally, uniformly and 
asymptotically. Compared to the linear counterpart, they 
accelerate the convergence and can provide for finite-time 
or fixed-time stability. 

I. Introduction 

This work is an extended version of [1]. In this 
paper the proof of the claims in [1] are presented in 
the appendices. This proof are omitted in the first work 
because of its length, instead, the space is used to 
discuss and clarify the results. 

A classical linear parametric model is given by 
y{t) := vF{t)9Q, where do € K" are the parameters, 
u{t) G M" is the regressor and y(t) € M is the measured 
signal. Output y(t) along with an estimate of the 
parameters 0(t) is used to build the output estimation 
error e(t) := u^{t)6(t) — y{t), which can be rewritten 
as e{t) = vF{t)x{t), where x{t) := 9{t) — do is the 
parameter estimation error. The aim is to use e to 
drive x to zero. Since Persistence of Excitation (PE) 
of u{t) is equivalent to the Uniform and Complete 
Observability of the associated linear dynamical system 
do = 0, y{t) = u^(t)9o [2], [3], it is a necessary 
condition to assure uniform and robust convergence of 
any algorithm. In particular, for the Linear Gradient 
Descent and the Recursive Least Square Methods PE 
is a necessary and sufficient condition for exponential 
convergence [2], [3]. 

In fact, the use of correction terms linear in e cannot 
provide for convergence faster than exponential. So, if 
accelerated convergence is desired, algorithms using 
nonlinear correcting terms are required. In recent years, 
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the use of homogeneous systems and homogeneous 
correction terms for control and observation purposes 
have been very successful in providing finite time 
and fixed time convergence [4], [5]. Eurthermore, 
homogeneous higher order sliding modes (HOSM) 
provide for discontinuous correction terms, which 
provide not only finite time convergence but also 
insensitivity to (matched, bounded) perturbations [6], 
[7], [8]. Homogeneity has been important for these 
results, since it provides useful properties: e.g. local 
asymptotic stability is equivalent to global finite time 
stability for systems of negative homogeneity degree 

[9], [6], [10]. 

Motivated by these results we propose in this 
work non-linear estimation algorithms, with non¬ 
linear correction terms in e, that lead to time-varying 
dynamical systems. The proposed schemes can be 
obtained as the negative gradient of an instantaneous 
convex non smooth function of the output estimation 
error. The resulting (error) systems are homogeneous 
(or homogeneous in the bi-limit) in the estimation error 
X [9], [11], but time-varying. Unlike time invariant 
homogeneous systems, time-varying systems do not 
possess such strong properties. Eor example, neither 
local asymptotic stability implies global asymptotic 
stability, nor negative homogeneity degree implies finite 
time convergence [11]. Some local asymptotic stability 
results for time-varying systems homogeneous in the 
state have been obtained, when the homogeneity degree 
is zero [12], and for positive homogeneity degree [11], 
using averaging techniques. 

We show in this paper, in straightforward a direct 
manner, that the proposed algorithms converge globally 
and asymptotically under the PE condition of the regres¬ 
sor. When n = 1, the algorithm is able to converge in 
finite-time and it is also able to estimate a time-varying 
parameter. Moreover, adding homogeneous terms of 
positive degree the estimation error converges in fixed¬ 
time, i.e. the convergence time is upper bounded by a 
constant independent of the initial estimation error. Eor 
n > 1 global, uniform and asymptotic stability can be 
assured and acceleration for large initial conditions can 
be obtained. 



II. Motivation and Problem Statement 

The estimation problem may be regarded as a mini¬ 
mization problem. In this context a common cornerstone 
is to set a convex cost function of the output estimation 
error. In this work we choose the following structure for 
the cost function 

Here the exponent p > 0 is a parameter to be chosen. 
The term —^J{d)p shows the direction in which the 
parameter estimated needs to change. With this idea in 
mind we propose the following algorithm 

'm ■■= 

oO 

= -\u^{t)e{i) - {t)e{t) - y{t))u{t). 

For the sake of readability let us define := 

|w|^’sign(w). With this convention the algorithm is 
rewritten as 

e{t) = - \u^{t)e{t) - y{t)\Pu{t). (1) 

We denote as composite algorithm an algorithm that 
results from adding vector fields of the form in Q to 
avail of the dynamics traits associated to each pi. This 
leads to 

u(t)- ( 2 ) 

To analyse the convergence of the algorithms to 9q the 
dynamics of the estimation error x{t) = 9{t) — 9 q 
is needed. Hence it is necessary to compute the time 
derivative of x{t). In the case of the single algorithm in 
Q this error dynamics is as follow 

x{t) = — \u^ {t)x{t)\Pu{t). (3) 

Considering f{t,x) = —\u^{t)x\Pu{t) and replacing 
the argument x by ex, with e > 0, we have the 
following relation f{t,ex) = e^“^e/(f,x). Taking this 
into consideration and following Definition 1 in [11], 
we can said that the system ([^ is homogeneous with 
homogeneity degree p — 1. When 0 < p < 1 the 
homogeneity degree is negative; for p = 1 it is zero; 
and for p > 1 the homogeneity degree is positive. Now, 
the convergence analysis for this class of algorithms 
reduces to establish the stability and the attractivity of 
the origin of a time-variant homogeneous system. 

Repeating the same analysis for the composite algo¬ 
rithm the error dynamics is 


Please, notice that this system is no longer homogeneous. 

A quick stability check with a quadratic Lyapunov 
function V{x) := yields to global uniform stabil¬ 

ity (GUS) of the ongin of ® and ffl. The time derivative 
of V is presented for both cases in order to show its 
negative semi-definitiveness: 

along (|^ V(t) = < 0, (5a) 

h 

along (|^ F(t) = — < 0. (5b) 

i=l 

Notice that the term u^{t)x{t) can vanish outside the 
set {t\u{t) = 0} U {t\x{t) = 0}. To assert the uniform 
asymptotic stability (UAS) of the origin of (|^ and 0 
u{t) needs to be PE. In [13] a convenient description 
of the Persistent Excitation is given. This description 
is expressed as a lower bound of the integral of the 
regressor and is found convenient because it fits well 
in the study of the estimation-error convergence. Eor 
convenience the definition is reproduced below 

Definition 1: Let u{t) : K+ —K" be a piecewise 
continuous function. It is said that u is of PE if there 
exist T > 0 and e > 0 such that 

t+T 

|u’^(s)u>|ds > e, 
for all in e K” with ||r(;|| = 1. 

A 

To conclude this section, we present some definitions of 
time-varying systems. 

Definition 2: Let a time-varying system be repre¬ 
sented by x{t) = f{x(t),t) where /(0,f) = 0 for all 
t. Let H be a connected open subset of K", such that 
0 € H. The point x = 0 is 

• Uniformly finite-time stable (UETS) if it is uni¬ 
formly stable and for any Xq S H there exist 
0 < T(xo) < -foo such that x{t,tQ,xo) = 0 for 
all t > to + T(xo). Also, if H = M" then x = 0 
is said to be globally uniformly finite-time stable 
(GUETS). 

• Uniformly fixed-time stable (UExTS) if it is GUETS 
and exist T < -foo such that T > T(xo) for every 
xo e M”. 

A 



III. Main Result 

In this section the stability of ([^ and Q is presented. 
Two cases are recognized; (i) when only one parameter 
needs to be estimated (scalar case) and (ii) when are 
more than one (vector case). This division is done 
because the results in the vector case do not reflect 
certain phenomena that occurred in the scalar one. 

The proof of the theorems can be found in the Appen¬ 
dices. 



A. Scalar Case 


The results in the scalar case are stronger than in the 
vector one due to the fact that the product {t)x(t) 
can only be zero if one or both of the variables are 
zero. Persistent excitation prevents u to stay in zero 
or to exhibit a growing dwelling time in zero. This 
implies that, for staying in zero, x needs to be zero. 
When the exponent p is chosen in the interval [0,1) the 
algorithm converges in uniform finite-time and a bound 
of the convergence time is given as an integer multiple 
of the persistent excitation period T. The next statement 
summarizes this discussion. 

Theorem 1: Let u{t) be a piece-wise continuous func¬ 
tion of t and of PE. Let 0 < p < 1 and n = 1, then the 
origin of © is globally uniformly finite-time stable. An 
upper bound of the convergence time is 


T(a;o) < 


(1 - p)TeP+i 


T 


( 6 ) 


where a;(fo) is the estimation error at the initial time; T 
and e are as in Definition and [•] denotes the ceiling 
function. 

A 

When p = 1 the linear case is obtained and its properties 
are well know [2], [3], for this is left out of the 
discussion. For p > 1 only global uniform asymptotic 
stability (GUAS) can be asserted but other interesting 
property arise. No matter how large the initial error is, 
the time t(c) need to reach a smaller level set V{x) = c 
is only function of c. This is referred to as escape from 
infinite in finite time uniformly in t. We gather these 
results in the following theorem. 

Theorem 2: Let u{t) be a piece-wise continuous func¬ 
tion of t and of PE, p > 1 and n = 1, then the 
origin of ([^ is globally uniformly asymptotically stable. 
Furthermore, the time needed to escape from infinity to 
a compact region V (x) < c is bounded by 


"(c) < 


2V(p- l)TeP+i 


C 2 


T (7) 


A 

The study of the convergence of the composite algorithm 
can be done using the previous results. For the scalar 
case the stability can be asserted via a Comparison 
Lemma for differential inequalities [14]. In short, the 
trajectories of the system Q are below the trajectories 
corresponding to each of the single algorithms with every 
isolated exponent pi. This can be seen from (|5^ which 
can be rewritten as 

h 

V{t) = - 

<-2^|u(f)|Pi+V^(f), i = ,h. 


An important case occurs when at least one exponent is 
in [0,1) and another one is greater than one. For this 
case the time needed for the algorithm to converge to 
9q is independent of the initial estimation error and the 
initial time if u{t) is of PE. This is summarized in the 
following theorem. 

Theorem 3: Let u{t) be a piece-wise continuous func¬ 
tion of t and of PE, also n = 1. Consider the system 
Q and let P := {pi,P 2 ,’'- ,Ph} be the set of the 
exponents. Denote pm the minimum element in P and 
Pm the maximum. Assume that 0 < Pm < 1 and 
Pm > 1- Let P^ be the subset of the exponent smaller 
than one and Pm the subset of the exponents greater than 
one. Define p G Pm as the exponent which maximizes 
{pi — l)eP'~^^ andp G P^ the exponent which maximizes 
(1 — Pi)ef‘+^, then an upper bound of the time needed 
to converge to zero is 


( 

1 

-f 

1 

1 

(1 - p)TeP+i 

2^(p- l)reP+i 


A 

Last but no least, a discontinuous algorithm capable 
of estimating one varying parameter is presented. This 
algorithm makes use of the regressor sign and the sign 
of the output estimation error 

§(t) = -L ■ sign(u(f)0(f) - y{t)) ■ sign(u(f)), (8) 

with L > 0. Now assume that the parameter variation is 
bounded, i.e. 6 G [—7,7], 7 > 0. Also it is considered 
that u{t) is of PE and cannot stay in zero for time 
intervals, but can cross it. The error dynamics induced 
by (|^ is 

x{t) G —L ■ sign(a;(f)) -f 0(f). (9) 

which is a differential inclusion [9]. By employing 
V{x) = ^\x\'^ as Lyapunov Function its derivative along 
the trajectories of (|^ is 

V{t) = —L|a;(f)| — 9x{t) 

= -{L + 0 ■ sign{x{t)))\x{t)\, 

which is negative if L > 7 and u{t) f 0. If ult) stay in 
zero the track of 9{t) is lost but can be recovered later 
if u is of PE. To guarantee exact tracking, u{t) cannot 
stay in zero. 

B. Vector Case 

When 9q is a vector the set where uf {t)x{f) is zero 
grows. This changes the general behaviour of our algo¬ 
rithms. Only GUAS can be guaranteed in general with 
persistent excitation. Also, the discontinuous algorithm 
which result of selecting p = 0 does not converge. In the 
latter case, a signal of PE can be constructed for which 
the output estimation error e(f) becomes zero in finite 















time but x{t) does not reach the origin. The stability 
properties of the origin of (|^ are summarized in the 
following theorem. 

Theorem 4: Let u{t) be a piece-wise continuous func¬ 
tion of t and of PE, uniformly bounded by um, then the 
origin of ([^ is globally uniformly asymptotically stable 
for any p > 0. 

A 

Although only GUAS can be asserted with PE in general, 
in the following section two classes of signals of PE 
are presented. One of them guarantees global uniform 
exponential stability (GUES) and the other GUETS; in 
both cases for 0 < p < 1. This means that the signals 
u{t) which can provide UETS are in a subset of those 
of PE. 

The composite algorithm still works in the vector case 
but nothing more than GUAS can be claimed. In contrast 
to the scalar case the stability of the composite vector al¬ 
gorithm cannot be obtained via the Comparison Lemma. 
The following theorem synthesizes this discussion 

Theorem 5: Let u{t) be a piece-wise continuous func¬ 
tion of t and of PE, uniformly bounded by um, then the 
origin of Q is globally uniformly asymptotically stable 
if Pm a Pm + 1 where p^ is the minimum exponent in 
the set {pi,p 2 , • • • ,ph} and pM is the maximum. 

A 

The extra condition regarding the exponents appears due 
to the way the proof was done and we think it is not 
intrinsic to the stability. 

In the scalar case the Comparison Lemma gives infor¬ 
mation about the relationship between the trajectories 
of the single algorithm but in this case nothing can be 
concluded. In the next section simulation examples are 
presented for the composite algorithm. In the Figure 
and Figure]^ it can be seen that the Lyapunov function is 
below of the corresponding one for the single algorithms 
but this does not need to be true in general. 

IV. Examples 

In this section examples are presented. In the scalar 
case, simulations of a parameter estimation process are 
shown, whereas for the vector case also the behaviour 
of the error is studied for a specific class of signals. 

A. Scalar case 

Numerical simulation is performed to illustrate the 
difference between the classic gradient algorithm and the 
family presented in this work. For this aim we choose: 

6*0 = 5 , 9{to) = 0, u{t) = 2 cos{2t). 

Four algorithms were simulated: one with pi = |, 
other with p 2 = 1, a third one with ps = § and a 
composite one with pi and pa. Figure and Figure 
show the behaviour of the Lyapunov function V{t) — 
f 10(f) — 00P- fn the first figure the initial value of 


V{t) and its decaying behaviour are presented for all 
the algorithms. As can be seen there is an ordering in 
the decay: Composite >p=|>p=l>p=|, and 
this is due to the initial value of U(fo) >> 1. When the 
estimation error becomes smaller the order changes as 
in Figure 0 

In Figurerathe behaviour is shown when V{t) < 0.014. 


V(t) 



Fig. 1. Behaviour of V{t) in the Scalar case (a). 

In this region the decay order is: Composite > p = 
|>p = l>p=|asis expected from the 
kind of convergence of every algorithm: u. fixed-time > 
u. finite-time > u. exponential > u. asymptotic. This 
order does not change for any future time. 

Now a estimation of a time-varying parameter is done 


V(t) 



Fig. 2. Behaviour of V{t) in the Scalar case (b). 

via The simulation parameters are: 

0(f) = cos(3f) — 4, 0(fo) = 0, u(t) = cos(f) -I- 1.5, 

L = 3.3. 

Since the regressor is non-zero for any time, an exact 
tracking of the parameter is achieved as shown in Figure 

[3] 

B. Vector case 

First simulation results are presented for the estima¬ 
tion process of three parameters with ([T]) and Q. Four 
algorithms were simulated: the first with pi = 3/4, the 
second with p 2 = 1, the third one with pa = 3/2 and the 
last one with pi and pa. The simulation results are shown 
in Figure]^ and Figure]^ to compare the behaviour of the 
algorithms, when different exponents are chosen. Later, 
the response of the algorithm for 0 < p < 1 is found for 
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e(t) 



Fig. 3. Time-varying parameter estimation process. 



Fig. 5. Simulation results vector case (b). 


piecewise constant signal and the class of convergence 
is established for it. 

For the simulation example the next conditions were 
selected: 

u{t) = [2 cos(2t) — cos(3f) 5 cos(5f)] ^ , 

00= [-3 ^/2 4]^, 

0(to)=[O 0 0]^. 

As in the last section only the plot of V(t) = — 

Oq)'^{ 9{t) — 9 q) is presented. In Figure W the decay of 
V{t) is shown at the beginning of the process and the 
decay order is preserved. Also in Figure the change 
of order that was found for the scalar case is obtained. 
However we do not see, in general, a change of order 
for arbitrary parameter values. 

Now the response of the system when u{t) is 


V(t) 



Fig. 4. Simulation results for vector case (a). 


piecewise constant is analyzed. Let us define a new 
variable as follows z{t) := u^{t)x{t), its dynamics is 
z{t) — u^{t)x{t) = —\\u{t)\\'^\z{t)\P since u{t) — 0 
almost everywhere, where || • || denotes the Euclidean 
norm. For simplicity we assume that the length of the 
intervals where u remains constant has a constant value 
r > 0. The solution of z(t) in the interval t G [ti,ti+T] 
with u{t) = /r, um > HfII > 0 and p > 0, p ^ 1 is 

z{t) = (| 2 (ti)|^“P - (1 - (i - D) ^ sign(z(ii)). 

For 0 < p < 1 the expression above holds if 
> (1—p)||u|p(f —fi) and z{t) = 0 otherwise. 
The solution of z{t) can be used to find the solution of 


x{t) by noticing that x{t) = —\z{t)\Pu. This is given 
in the forthcoming expression 


x{t) = - sign(p^a;(ti))x 

Restricting the analysis for 0 < p < 1 it is clear that, if 
T is large enough then x{ti +r) is orthogonal to 

x(ti + r) = x{ti) - — 


Fix T, if ||x(fi)||^“^' < (1 — p)u^^ T then x{ti + r) 
becomes orthogonal to /i. This means that there exists 
a ball centered in zero for which r is always large 
enough to make x orthogonal to any p. If the sequence 
U = is chosen to fulfil the notion of persistent 

excitation for discrete systems in [15], then the origin 
is GUAS by Theorem and x{t) can reach any ball 
centered in zero in finite time, i.e. always reach the ball 
in where t guarantees that x becomes orthogonal to pi, 
and this yields a discrete system which is described by 
the following difference equation 

Xk+l = fin ~ T] 

IlFfeIr > 


Only assuming PE of the sequence [/, GUES of the 
origin can be concluded making an analogue analysis to 
the o ne shown in the proof for Theorem 4 in Appendix 
II-D| but taking 14.(x) = x\xk, AV = Vk+r - Vk{x) 


instead of V(x) and V{t) respectively. Now take n 
mutual orthogonal vectors {ui,U 2 , • • • ,fn}. 'vjvj = 0, 
i ^ j, and excite the system with them, then Xk+n = 0 






• Xk 



E 



Xk = 0. 


A sequence constructed with these vectors is of PE and 
also makes the system GUETS. This shows that persis¬ 
tent excitation cannot guarantee finite-time convergence 
in the vector case but it does not forbid it. 





































V. Conclusions 

In this work a parameter estimation technique is pre¬ 
sented. With the proposed algorithms we obtained hnite- 
time and hxed-time convergence to the true parameters. 
However this properties cannot be guaranteed in general. 
A deep study of the signals that can assert such important 
properties is still needed. 

Even though the algorithms were selected to make the 
error dynamics homogeneous in the state, this does not 
help in the analysis. However, the homogeneous non- 
linearities can enhance the robustness properties of our 
algorithms w.r.t additive perturbations in comparison 
with classical approaches. 
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Appendix I 

Persistent Excitation 

To prove uniform asymptotic stability rather than 
uniform stability the persistent excitation in the regressor 
is needed. For using the property adequately the next 
proposition is developed. 

Proposition 1: If u{t) is of PE, then the following 
inequality holds for T, e and w as in Definition and 

p > 0 

pt-\-T 

Proof: Applying the Holder inequality to f{f) = 
\u^{t)w\ and g{t) = 1 on the interval [t, t+T] and using 
and as Holder conjugates we obtain 


ft+T / j-t+T \ 

y |u^(s)ui|’’'''^ds > Ij |u^(s)u;|dsj 

using Dehnition leads to 

/ t-\-T rpp-\-\ 

|u^(s)w|P+Ms > 


P+1 


.p+1 _ J^^p+l 


This derivation from the PE is done in the aim of easily 
present the proof of the theorems. 

Appendix H 

Proof of Theorems[I]to[5] 

A. Proof of Theorem 

For n = 1 the term \u{t)'^x{f)\ can be rewritten as 
|u(f)| • \x(t)\ — s/2\u{t)\V^ (t) replacing this in ( |5al l 
yields 

V{t) = -2^|M(t)|P+iv'7^(f). 

Solving the differential equation we have 


V{t)= (v5^{to)-2"-+il-p) 


|w(s)|^+Ms 


*0 


This solution is valid for p > 


( 10 ) 

1- 


0 and p 

For 0 < p < 1 the solution exists if V^ifo) > 
2^ (1 —p) |M(s)|^+^ds after that V{t) = 0. The PE 

guarantees that there exists ti > tg when the inequality 
no longer holds. To estimate this time it is sufficient to 
find an integer k such that the integral of from 



to to to + kT is greater than 2 (to)- 

Using Proposition [T] 


f'tQ-\-kT 


1—p 

|m(s)|p+Ms > kTeP+^ > 


1-p 


Solving for k and taking the least integer that fulfills the 
inequality yields 


k = 


(l-p)reP+i 


V^ito) 


( 11 ) 


then the time that guarantees that y(t) reaches zero is 
ti = k ■ T + Iq. Notice that U(to) can be replaced by 
i|a:(to)P in ([TT]l to obtain (|^. 


B. Poof of Theorem 

Consider again the solution for y(t) in ( [TOl l. Since 
p > 1, Equation ( [TOl i can be rewritten as 

V{t) = --^ 

(u'2’’(fo) + 2'’2 '(p- 1)//^ |u(s)|P+idsj ” ' 

As u{t) is of PE the denominator grows unbounded 
making V{f) —>■ 0 as f —>■ oo. Now an estimate of the 
time needed for V (t) to decreases from V (to) to a value 
equal or smaller than c is calculated. Substituting V{f) 
for c in ( |5a| i and evaluating the integral from to to ti it 
is clear that the value of the integral needs to be larger 
enough to satisfy the next inequality 
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Eixing ti — to as an integer multiple of T, i.e. ti — 
to = kT, from Proposition we know that the integral 
is greater or equal to kTe^'^^. Eorcing the RHS of the 
last inequality to be less than kTe^'^^ we get 
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Now solving for k and taking the smallest integer that 
fulfil the inequality yields 
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Taking the limit when U(to) —oo the bound Q is 
found 
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C. Proof of Theorem 

Prom ( |5b| i it follows that Vit) < —\uF{f)x{f)\P'^^ = 
— 2 ^ 2 —|u(t)|U^ 2 —(f), which define a different differ¬ 
ential inequality for each term in the sum. Prom the 
Comparison Lemma we known that the solution of V (t) 
is for below of each solution of Vpi(t), where Vp.{t) = 

—2^^ |u(f)|Uj ^ (f); Vp.{t) take the form of ( [T0| ). Let 
Pm and be as in the theorem. Por each pi G Pm 
we can assert that Vp^{f) can escape from infinity to a 
compact set in finite time and so V(t). Pix the level set 
as V(x) = 1 and estimate the time needed to reach it 
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The smallest time that the algorithm can guarantee is 
obtained when {pi — l)eP’+^ is maximum. Take that 
quantity as the estimate. Now, with pj G Pm, we can 
estimate the time needed for each Vp- (<) to converge 
from the level set V (x) = 1 to zero 


T(a;o) < 
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T, Pj G Pm,xo\V{xo) = 1. 


Again, the smallest time that the algorithm can guarantee 
is when (1 — pj)e^ 2 +^ is maximized. Then the time 
needed by the algorithm to converge is, at most, the sum 
of the two estimates. 


D. Proof of Theorem 

Pollowing the idea in [16] a lower bound of the 
integral of \u^{t)x{t)\P~^^ is needed. Take the term 
lu^{t)x{ti)\ and add a zero in the form vf'{t)x{t) — 
u^{f)x{f) inside the absolute value, by means of the 
triangle inequality the next partition is obtained 

\u^{t)x{t)\ > \u^{t)xiti)\ - \u^{t){x{ti) - xit))\. 

Rising both sides to p +1 and using the Jensen inequality 
after that, the inequality yields 

\u^it)x{t)r^>^\u^it)xih)\p+^ 

- - x{t))\P+\ 

Integrating both sides from ti to ti+T and using Propo¬ 
sition [T] to bound the first term in the RHS becomes 

r-tj+T T^^P+l 

- \u^it){xih) -xitW+^dt. (12) 

Jtl 

To bound the magnitude of the second term in the RHS 
of the last inequality, the next procedure may be used 






























assuming u(t) is uniformly bounded, i.e. ||w(t)|| < um'- 
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The norm of x in this case is {t)x{t)\P\\u{t)\\ which 
is less than \u^{t)x{tM'‘’uM- Holder inequality can be 
used with /(t) = \u^(t)x{t)\P and g{t) = 1 in the 
interval + T] with and p + 1 as Holder 

conjugates, then 
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Using this in ( [T3] l and then the result in one gets 
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This can be rewritten as 


PgP+i 

2P 


Jtl 


By algebraic manipulation and defining z := 
{t)x{t)\P~^^dt the notation can be simplified; 

z+TVi"+v> ^iix(u)r+'- 

Notice that the polynomial P(z) := z + is 

a strict monotonically increasing function for z > 0 and 
its inverse exists. By denoting this as P~^{z) and using 
in the inequality above, one gets 

fti+T \ 

z=y^ \u^ > p~^ . 

Recalling V(t) from ( [Sa] ) and integrating it for the same 
time interval it is obvious that 


/ 

Since P“^() is also a strict monotonically increasing 
function and from Theorem 5 in [13] the global uniform 
asymptotic stability of the origin of 0 is stablished for 
any p > 0. 
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E. Proof of Theorem 

In ( |5b| ) several terms of the form \w^{f)x{f)f^'^^ 
appear; an analysis of each term separately is done 
before study the full derivative of U. By taking the term 
\v?'{t)x{ti)\ and adding a zero in the form vf {t)x{f) — 
u^{f)x{f) inside the absolute sign and then applying the 
triangle inequality the next partition is found 

\u^{t)x{ti)\ < \u^{t){x{ti) - x{t))\ + \u^{t)x(t)\, 


rising to pi + 1 and using the Jensen Inequality 

it)x{ti)f'+^ <\u^it){x{ti) - 

+ \u^{t)x{t)\P'+\ 

Solving for {t)x{t)\P'^^ and integrating between 
[ti, ti + T] 

/ |«^(s)a;(s)|P>+Ms > — / |u^(s)3;(ti)|P'+ ds 

J t\ ^ ^ J t\ 

Pt-l+'l' ^ 

— / |ri^(s)(3;(ti) — a;(s))|P»”*"^ds. 
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A lower bound of the first term in the RHS of the in¬ 
equality can be found using Proposition [T] The following 
inequality is obtained 

rti+T 7^fPi + l 
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(14) 

A upper bound of the magnitude of the second term 
is quite more difficult. Assuming u{t) is uniformly 
bounded by um, i-C- ||ii-(i)|| < um, the integral can be 
bounded as 
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Now an analysis of the integral of ||i(f)|l is required. 
First the norm of x is estimated 

h h 

||i|| = I ||u(i)|| < UMy^fu^{t)x{t)fJ. 

i=i 

Integrating and using the Jensen inequality for the 
function 

ftl+T NPi + 1 
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The restriction pM < Pm + 1 guarantees that > 1 

Pj 

for every (*, j) and then are Holder 

conjugates. Now the Holder inequality is used with 
f{t) = lu^(t)x(t)\P^ and g{t) = 1 in the interval 
[ti,ti +T] and Holder conjugates as described before 

. p^+l 

|ii^(s)a:(s)|^J dsj 

(17) 

Using ( [T5] | in applying ( [T7| l in the sum, one has 

„2^P. + 1)^P.^P,+2 j- ^ + 

||i(s)||ds). 
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From ( [T4l i, ( fTSj l and ( [T8] l we obtain 
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Defining Zj := |M^(s)a;(s)|7'‘+^ds the last in¬ 

equality can be rewritten as 
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The polynomial in the LHS is a strict monotonically 
increasing function and zero when = 0, then its 
inverse exists and also is a strict monotonically increas¬ 
ing function. By denoting the polynomial as Pi and its 
inverse as P~^ one gets 

r^i+T \ 

> 77- —||a;(ti)||P‘+lj . 

(19) 

Now integrating from fi to fi -f T and using ( [T^ 
the next inequality for V is obtained 
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This satisfies the assumptions of the Theorem 5 in [13] 
and guarantees the global uniform asymptotic stability 
of the origin of (|^. 







